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Abstract 

The relationship between measured transverse energy, total charge recovered in the detector, and 
size of the emitting system is investigated. Using only very simple assumptions, we are able to 
reproduce the observed binomial emission probabilities and their dependences on the transverse 
energy. 
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During the last decade, evidence has been mounting that nuclear matter undergoes a 
phase transition in the nuclear fragmentation process. From general consideration regarding 
the elementary nucleon-nucleon interaction (repulsive at short and attractive at intermediate 
distances), we expect the nuclear phase-diagram to show a Van-der-Waals “liquid-gas” phase 
transition of hrst order, terminating in a second-order transition at the critical point. 

Recent observations point to evidence of hrst and second order transitions. In experi¬ 
ments studying Au-Au collisions conducted at the GSI, a measurement of the temperature as 
a function of excitation energy found possible evidence for a two-phase coexistence regime [|l| , 
not unlike the scenarios predicted by statistical multifragmentation models with excluded 
volume PQ. Other experiments conducted at the Bevalac focussed on the extraction of 
critical exponents from (almost) completely reconstructed Au-fragmentation events on C- 
targets, studying the dependence of the second moment of the charge distribution and size 
of the largest fragment as a function of the total charged particle multiplicity H- It was 
shown that these data are consistent with the second-order phase transition predicted 
by the nuclear percolation model P HT^] . 

If one wants to gain a fundamental understanding of the fragmentation process that goes 
beyond simple equilibrium model descriptions of the phenomena, then a proper description 
of the origin and time evolution of fluctuations is essential ||13|-|I^ , in particular if one wants 
to understand why particular molecular dynamics codes produce fragments (or not!), and 
what their connections to the fundamental processes of nuclear fragmentation are ||T7|-p7|. 

In this light, the recent hndings of Moretto et al. are all the more surprising |28|-|3^. 
This group found that the probability Pn of emitting n intermediate mass fragments (IMFs) 
follows a binomial distribution 

Tfl ^ 

Pn{m,p) = • p-(l-p)—- (1) 

nl [m — n)l 

The parameters m and p are related to the average and variance of the distribution. 

C50 

(n) = '^nPn {m, p) = m ■ p (2) 

n=0 
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( 3 ) 


(T^ = ^ (n — {n)Y Pn{'m,p) = rnp{l — p) 

n=0 

This result suggest that one may interpret the parameter p as the elementary probability 
for the emission of one fragment and the parameter m as the total number of tries. This 
would indicate that the problem of multi-fragment emission is reducible to that of multiple 
one-fragment emission. The claim for reducibility and its interpretation as the consequence 
of a simple barrier penetration phenomenon was further strengthened by the observation 
that ln(p“^) has a linear dependence on 1 /a/^, where Et is the total transverse energy, 
Et = J2i Oi. Finally, the same scaling was found for different beam energies and 

different projectile-target combinations. 

Other authors have criticized the above work, pointing our that there are different emis¬ 
sion probabilities for different size IMFs, that there are problems in the transformation 
between the total transverse energy and a true thermal energy |^, and focussing on au¬ 
tocorrelations between the number of IMFs and the transverse energy [Q. (See also the 


replies to these criticisms in refs. 


In the present note we add to this discussion by showing how binomial distributions 
arise naturally from hnite size effects. In particular, we focus on the dependence of the 
experimentally recovered charge as a function of the measured transverse energy. We then 
demonstrate why the dependence of the binomial parameter p on the total transverse energy 
arises. 

We begin our study by generating power-law distributed random fragmentation events. 
This is accomplished by determining the charge of individual fragments with a probability 
distribution proportional to Z~'^, where Z is the fragment charge, and r is the power-law 
exponent. For dehniteness, we wish to generate events with exactly Z^yg charges. If an 
event has less than Zgys charges, we add another fragment; if it has more than Zsys charges, 
we throw it out. For an inhnite system, we would expect the multiplicity distributions for 
individual fragments of a given Z to follow a Poisson distribution. 


^ , A” exp(—A) , , , , , 9 

Q„(A) =-j- with : A = (n) = 


n\ 


( 4 ) 
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And since the combined probability distribution of two Poisson-distributed variables is again 
a Poissonian, 

n 

Qn{^l) ® Qn{^2) = X/ ■ Qn-i{^2) = Qn{^l + A 2 ) , (5) 

i=0 

we would expect that the multiplicity distribution of the total number of intermediate mass 
fragments (IMFs) is also Poissonian. 

The individual probability distributions for IMFs, however, cannot be exactly Poissonian, 
because the tails of the distributions are cut off due to the hnite size of the emitting system. 
Thus the probability distributions in our simulation are closer to a binomial distribution 
with rather large values of m and small values of p. (When m ^ oo, p ^ 0 such that 
mp =const., we obtain a Poissonian as the limit of a binomial distribution.) Typical values 
of p we hnd for the probability distributions of our individual fragments are < 3 ■ 10“^ for a 
system of 100 total charges. These small values of p imply that the probability distributions 
are very close to a Poisson distribution. 

We now ask what the combined probability distribution for fragments charges in the 
interval 3 to /c, /c = 4, 5,..., 30 is. (If we use k = 20, this corresponds to the usual dehnition 
of IMFs.) We hnd numerically that to very good approximation this distribution is again 
a binomial distribution, for all values of k. The binomial parameters pk and rrik have 
a monotonical behavior as a function of k: pk rises monotonically until it saturates at 
k = and ruk falls monotonically. This results directly from the mathematical fact that 
the mean of the combined probability distribution is the sum of the mean values of the 
individual distributions, but the variance is always smaller than the sum of the variances 
for the folding of binomial distributions. 

In hgure 1, we show the behavior of the parameters pk and of the combined probability 
distributions as a function of k, the upper limit charge for the folding procedure, for different 
values of Zgys- For each value of Zgys, we generated 10"^ events. This hgure already contains 
the essential key to understanding the patterns observed by Moretto and collaborators. We 
can clearly see that as we include more and more fragments in the dehnition of IMFs the 
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extracted values of the binomial parameter pk increase, and those for ruk decrease. We can 
also see that ps (= Pk{k=3)) decreases as we increase the total charge of the fragmenting 
system. This is expected: the larger the total available charge, the closer the probability 
distributions of individual fragments (in this case Z = 3) will be to the Poissonian limit. 
It is essential to note that the values of pk for each system size saturate a.t k = Zgyg. This 
is obvious, because we cannot have IMFs larger than the total charge available. But this 
obvious fact has an interesting consequence: The smaller the system size, the fewer the 
terms that can contribute to the construction of the asymptotic value of pk, and the lower 
the asymptotic value of pk- 

This fact, combined with the dependence of Z^ys on the transverse energy, already is the 
explanation for the scaling observed by Moretto et al. 

What is the dependence of Zsys on the transverse energy, Et, in the experiments of 
Moretto et all This is shown in £g. 2 for the reaction Kr+Au at 55 AMeV p6|,|T7|. The 
hlled plot symbols show the mean Z^ys for each value of the Et, and the error bars give the 
width (standard deviation) of the distribution. The dominant feature of this figure is the 
linear rise of the mean value of Z^y^. with Et, 


{Zsys{Et)) ~ 2 + 0.092 T/t/MeV 


( 6 ) 


for values of Et less than 0.7 GeV, and the saturation of Zgys for larger values. 

The width of the Zsys(i77’)-distribution is significant, on the order of 10 units of charge. 
If we wish to construct the probability distributions of intermediate mass fragments by using 
our knowledge of the dependence of the binomial parameters p and m on the system size 
(fig. 1), and the dependence of the system size on transverse energy, we have to integrate 
over the experimentally measured width of the Zsys(i7r)-distribution. The resulting values 
of p for the integrated distributions are shown as a function of Et in £g. 3. We have run 
three different calculations, using three different values of the exponent r in our fragment 
production probability distributions, p{Z) oc Z~'^. In this figure, we display the results in 
the same way that Moretto et al have done. One can already see that there is qualitative 
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agreement with the tendencies observed by the Moretto-group: for all values of r, we observe 
an approximately linear rise of lnp“^ with 1 /a/E^. This, however, is not the consequence of 
some kind of thermal scaling. Instead, it is purely a consequence of the variation of the size 
of the emitting system as a function of the transverse energy, and with it a change in the 
effective parameter p in the binomial probability distribution. 

The projectile and target masses only enter into our consideration as upper cutoffs for the 
possible maximum values for Zgys, and with it the upper values of Et. The experimentally 
found functional dependence of Zsys on Et that we show in £g. 2 for the system Kr+Au is 
basically the same for all target-projectile combinations and beam energies below approxi¬ 
mately 100 A MeV (for higher beam energies, we have sizeable radial flow contributions to 
the transverse energy). Transverse energy basically measures impact parameter, that is to 
say size of the emitting system; and for any experiment measuring inclusive fragment dis¬ 
tributions the results will be very much similar to the ones displayed in £g. 2. This explains 
the universal scaling observed by Moretto et al. without the need for invoking some deeper 
reason for this apparent universality; a plot of lnp“^ vs. l/^/E't is dominated by effects of 
the variation of the size of the emitting system. We should point out that our hndings are 
not dependent on the fact that the intermediate mass fragments carry transverse energy 
themselves. The only correlations entering our analysis are the experimental ones between 
transverse energy and system size. 

We can obtain more-or-less complete agreement with the experimental data, if we allow 
the power-law parameter r for the fragment mass distribution to vary with impact parameter 
and with beam energy. This variation is a well-documented experimental fact p7H3^; the 
experimentally observed value of r increases with impact parameter and therefore falls with 
transverse energy. If we assume 


T{Et) = 3.5 - Ei/(0.5 GeV) , (7) 

then we get the result displayed in fig. 4. Our calculations are represented by the plot 
symbols. The error bars are statistical and computed on the basis of 2 x 10^ events for 
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each point. The solid line is a ht to the experimental results of Moretto et al. As one 
can see, there is very good agreement. Assuming other functional dependences of r of 
may even yield better results. This agreement, however, is not quite as relevant as the 
main message we wish to impress on the reader: The universal scaling of lnp“^ vs. l/^/El 
is almost exclusively due to the hnite size of the system emitting the fragments and the 
dependence of the measured value of the transverse energy on that size. 

Even though the main message of the present note is that the Et-dependence of the 
extracted binomial parameters of the fragment multiplicity distributions can be explained 
rather straightforwardly, we do not wish to convey the message that there is no interesting 
information that one can extract from this type of analysis. For instance, the effects of 
varying system size could be eliminated with utilization of completely reconstructed frag¬ 
mentation events. For these types of events, percolation models predict a transition between 
sub- and super-Poissonian fluctuations near the percolation threshold . This type of be¬ 
havior is not expected in a sequential model. Once the kind of correlations discussed by us 
above are removed, then this type of fluctuations analysis should yield insightful information 
about the character of the nuclear fragmentation phase transition. 
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FIGURES 
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FIG. 1. Dependence of the binomial fit parameters pk and nik of the probability distributions 
for intermediate mass fragments from charge Z = 3 to Z = /c on the upper summation limit k, for 
different total charges of the fragmenting system, Zgyg. 
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FIG. 2. Dependence of the total charges of detected in a fragmentation event, Zgys, on the total 
transverse energy, Et, detected in the experiment 55 ^MeV Kr+Au |^]. The error bars indicate 
the width (standard deviation) of the distribution on an event-by-event basis. 
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FIG. 4. Dependence of the binomial parameter p of the IMF distribution on the transverse 
energy, assuming that the effective power r of the fragment probability increases linearly with 
transverse energy. Plot symbols with error bars represent our calculations, the solid line is a fit to 
the experimental data of Moretto and collaborators. 
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